Nonperturbative effect of attractive forces in viscous liquids 
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We study the role of the attractive intermolecular forces in the viscous regime of a simple glass- 
forming liquid by using computer simulations. To this end, we compare the structure and the 
dynamics of a standard Lennard- Jones glass-forming liquid model with and without the attractive 
tail of the interaction potentials. The viscous slowing down of the two systems are found to be 
quantitatively and qualitatively different over a broad density range, whereas the static pair corre- 
lations remain close. The common assumption that the behaviour of dense nonassociated liquids is 
determined by the short-ranged repulsive part of the intermolecular potentials dramatically breaks 
down for the relaxation in the viscous liquid regime. 
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Differentiating the respective roles of repulsive and at- 
tractive intermolecular forces in the properties of flu- 
ids and liquids has a long history in statistical mechan- 
ics [lj. The so-called 'van der Waals picture of liq- 
uids' H, Hi > *- e - the predominance of the short-ranged 
repulsive part of the intermolecular potentials in deter- 
mining the structure of dense nonassociated liquids, has 
proved very fruitful for predicting the pair correlation 
functions and the thermodynamics. Although not as 
thoroughly tested Jl, d, @, 0] , it has been useful for the 
dynamics as well [8[. 

More recently, this picture, in which the attractive part 
of the interactions is treated as a mere cohesive back- 
ground amenable to perturbative treatment, has been 
transposed to the viscous (supercooled) liquid regime. A 
number of approaches either suggest or take for granted 
that the structure and the dynamics of viscous glass- 
forming liquids are controlled by the short-ranged repul- 
sive forces. Among them are studies based on the mode- 
coupling theory of glasses [§] , the self-consistent phonon 
theory [id ], model potential energy landscapes [ll|, as 
well as recent work focusing on the correlations between 
pressure and energy fluctuations 12 j_ T^] and on the den- 
sity scaling of the relaxation time [l4| Predominance of 
the short-ranged repulsive forces is also shared by the 
'jamming scenario', which postulates that the physics 
of glasses and glass-forming liquids is controlled by a 
zero-temperature critical point characteristic of the jam- 
ming 15, 1^1 or glass 17 1 transition of spheres with finite- 
ranged repulsive interactions, with the longer-ranged at- 
traction considered as a perturbation. 

However, and quite surprisingly, the central hypothe- 
sis underpinning this large set of ideas has never been 
directly studied. It therefore appears timely to assess 
the role of the attractive intermolecular forces in the vis- 
cous liquid regime. Our central conclusion is that in the 
regime that is probed experimentally, attractive forces 



have little effect on the structure of the liquid, but affect 
their dynamics in a strong, likely nonperturbative, way. 

In this work, we compare by means of Molecular 
Dynamics simulation the structure and the dynamics 
of a standard model of glass-forming liquid, the Kob- 
Andersen 80:20 binary Lennard- Jones mixture (l8j and 
its reduction to the purely repulsive part of the pair 
potentials proposed by Weeks, Chandler and Andersen 
(WCA). In the following, the former is denoted by 'LJ' 
and the latter by 'WCA'. The interatomic pair potential 
between species a and (3, with a, (3 = A, B is given in the 
two systems by 



v a p{r) = 4e Q/3 



r < r 



and is zero otherwise, where r a/3 is equal to the position 
of the minimum of v a p (r) for the WCA potential and to a 
conventional cutoff of 2.ba a p for the standard LJ model; 
Cap is a constant such that v a p(r a n) = 0. The simula- 
tions are performed in the NVE ensemble (after equili- 
bration at a chosen temperature) with N — 900 — 1300 
particles (depending on the density) and with periodic 
boundary conditions. A broad range of density has been 
considered with p from 1.1 to 1.8. Lengths, tempera- 
tures and times are given in units of <jaa, ^aa/^b, and 
{ma\ A / 4S>eAA) 1 ^ 2 respectively. In line with the WCA 
theory, the two liquid models are compared at the same 
(p, T) state points. Their pressure then differs, with the 
attractive interaction roughly providing a negative back- 
ground term, as we have directly checked by studying the 
equations of states, P — P(p, T). 

We first consider the static structure of the liquids as 
characterized by the pair correlation functions g a /3( r )- 
The pair correlation of the density fluctuations, g(r) = 
x aXjig a ji(r) with x a the concentration of species a, 
is displayed in Fig. [1] for two temperatures at the typical 
liquid density of 1.2. As anticipated 0], the attractive 
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forces play virtually no role in the high-T liquid. The 
same is found at very high density (p > 1.6, not shown 
here) for all T's. We observe that the contribution of at- 
traction remains very small over the whole range of (p, T) 
under study. Inspection of the partial correlation func- 
tions g a p{r) shows the same feature. A detectable effect 
is observed on the self-correlation of the minority com- 
ponent, gBB{r), but since xb = 0.2, it has little impact 
on g(r). 

Turning now to the dynamics, we present the data for 
the time dependence of the self intermediate scattering 
function 



F s (q,t) = 




iq-(rj(*)-rj(0)) 



(1) 



with qa aa — 7.2, which corresponds to the position of 
the the peak of the total static structure factor at p = 1.2. 
We see that in the high-T liquid (Fig.fTJa), the agreement 
between WCA and LJ is not perfect but is quite good (a 
factor of 2 or less difference in the relaxation time for 
T > 1, p = 1.2). However, the difference between the 
two systems rapidly increases as one enters the viscous 
regime and becomes enormous at the lowest accessible 
temperatures. The relaxation time r a of the WCA model 
is then more than 3 orders of magnitude faster than that 
of the LJ model (see Figs.[Tb and[Tb). 

This is confirmed by other measures of the slowing 
down of relaxation. For instance, the empirically de- 
termined mode-coupling singularity temperature T c , ob- 
tained by an algebraic fit of the temperature dependence 
of r a , is roughly divided by two when removing the at- 
tractive forces. We point out that, although not appre- 
ciated before, similar conclusions can be drawn by com- 
paring the results already published in the literature for 
different Lennard- Jones mixtures 18, [3] and for their 
WCA truncations [U[2lj]. 

We find that it is only by going to very high densities 
that the difference shrinks and the relaxation times of 
the two systems become more comparable. Yet, even at 
a density as high as 1.6, i.e. more than 30% above the 
commonly used density of 1.2, there is almost an order of 
magnitude difference between the relaxation times of the 
two systems, and at p = 1.8, relaxation times still differ 
by more than 20 % at low temperatures, despite the fact 
that the pair correlation functions virtually coincide at 
these densities. These observations prove that the dy- 
namics is not primarily determined by the static pair 
correlation functions, for instance raising some doubts 
on the ability of the mode-coupling theory to properly 
describe the phenomenon. 

The above results unambiguously show that the relax- 
ation of the LJ is considerably slowed down by the pres- 
ence of the attractive forces as the liquid is cooled down 
in the viscous regime and approaches the glass transition. 
One may however wonder whether this is a mere quanti- 
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FIG. 1: Comparison between the LJ and its WCA description 
at the typical liquid density p = 1.2 for several temperatures, 
(a) Static pair correlation functions g(r) as a function of r; (b) 
Time dependence of the self-intermediate scattering function 
for the majority component A; (c) Arrhenius plot of the relax- 
ation time T a defined from F 3 (q,t = r a ) = 1/e (the lowest T 
is 0.43 for LJ and 0.26 for WCA), with high-T Arrhenius fits 
shown as dashed lines. Despite similar pair structure, both 
systems display dramatically different dynamics. 



tative effect that can be accomodated by introducing an 
effective energy scale that takes into account the mean 
influence of the attractive part of the potentials by renor- 
malizing the temperature scale 



14j,|2J]. To test this hy- 



pothesis, as well as to represent all of our data for differ- 
ent densities on the same graph, we have fitted the tem- 
perature dependence of the relaxation time in the high-T 
regime to an Arrhenius formula: r Q ~ Too exp(_E 00 (p)/T). 
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This allows us to collapse the data at high T with a good 
accuracy and to extract an effective activation energy 
scale E 00 (p). The latter can then be used to compare 
the relaxation data in the presence and in the absence of 
attraction on a renormalized temperature scale T/E oa (p) 
and therefore to test the above hypothesis. No emphasis 
is put on the physical meaning of this Arrhenius fit, which 
we take as a convenient and nonsingular representation 
of the high-T data. 

The results are shown in Figs. [2^, Eb where we plot 
the logarithm of r Q for LJ and WCA for all densities be- 
tween 1.1 and 1.8 as a function of the inverse of the scaled 
temperature, E oc (p)/T. By construction all curves co- 
incide at high (scaled) temperature above some 'onset', 
T/Eoo(p) ~ 0.3, at which the viscous regime roughly 
starts and departure from simple Arrhenius fit becomes 
significant. Below this onset temperature, we find that 
all LJ data essentially collapse onto a master curve (with 
a small deviation seen for the lowest density of 1.1), as 
roughly do the WCA data for the three highest densities 
(p = 1.4,1.6, 1.8). The coincidence between LJ and WCA 
rescaled data is only fair at those densities, and the curves 
clearly diverge as one lowers the density to reach values 
more typical of regular supercooled liquids, i.e. p = 1.2 
(compare Figs. [2^ and [2}}). The isochoric 'fragility' of 
the WCA model is strongly density-dependent, which is 
reminiscent of the behaviour found in dense fluids of har- 
monic repulsive spheres 17| . This is clearly at variance 



with the almost constant isochoric fragility of the LJ. 

For completeness, we display in Fig. [2fc the density 
dependence of E x (p) for the 2 systems in a log-log rep- 
resentation. For LJ, it roughly goes as p 5 , in agreement 
with previous work finding relaxation data collapse with 
the scaling variable p 1 /T with 7 ~ 5 fl4j |. 

The viscous slowdowns of the LJ and WCA models 
are therefore not only quantitatively different at a given 
density, they are also qualitatively distinct. The density 
scaling of the relaxation that is empirically found in real 
glassforming liquids and polymers 23j, |2J, |25|, 12611 . and 
in the LJ model as well (see Fig. [2^ and Refs. [l4 . [23}) 
is strongly violated when attrative forces are truncated. 
These findings show that, contrary to expectations, the 
attractive components of the pair potentials play a cru- 
cial role in the viscous liquid regime when approaching 
the glass transition. A purely repulsive WCA system of 
course displays a slowing down of relaxation that should 
end up in glass formation at low enough T, but some 
of the characteristics of this slowing down, including the 
absence of density scaling of the relaxation time, are at 
odds with the behavior of the full LJ model it is supposed 
to describe, and of real glass-formers. 

Finally, we map out in Fig. [3] the various regimes stud- 
ied here in a (p, T) phase diagram and discuss the rel- 
evance of our findings to real glass-formers. On top of 
the thermodynamic transition lines, we have plotted the 
empirically determined mode-coupling line, T c (p), as an 
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FIG. 2: Rescaling of the relaxation data for the LJ and WCA 
models over a wide range of densities and temperatures. We 
use the activation energy scale Eoo (p) obtained by fitting the 
high-T data to an Arrhenius formula, (a) Arrhenius plot of 
the relaxation time for a scaled temperature T/Eoa(p) for LJ; 
(b) Same plot for WCA. Note the large change of fragility 
with density for p < 1.4, not seen in the LJ data, (c) Density 
dependence of the activation energy scale for the two models, 
with LJ data fitted with a power law. 



indication of the trend for the (isochronic) glass tran- 
sition line in the diagram. One can schematically dis- 
tinguish three regions. Region (I), inside the gas-liquid 
coexistence curve (or spinodal [28|), can only be accessed 
by removing the attractive part of the potentials. This is 
the region that could be controlled by a zero-temperature 
jamming 15, li| or glass 17| critical point. Region (II) 
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FIG. 3: Phase diagram of the LJ model. The melting line 
(shown as a dashed line) is for the one-component model [27T ]. 
We distinguish 3 regions. (I) is inside the liquid-gas spinodal 
and is only accessible for WCA: see the zoom in the inset. 
The extrapolated curve shown as a dotted line ends in a pu- 
tative zero-temperature singularity [TBI [l7j. (II) is the range 
corresponding to the experimentally accessible liquid regime, 
where large differences are found in the dynamics of the two 
models. (Ill) is the high-T and high-p region only reachable 
in computer simulations of model systems, in which the role 
of the attractive forces becomes small or negligible. 



is the experimentally accessible range of glass-forming 
liquids, for which, typically, the pressure can be varied 
from 1 bar to 10 kbars with an associated density varia- 
tion of 20 - 25%. Finally, region (III) is a high-T, high-p 
regime that is only reachable in computer simulations of 
model systems. It is only in this region that the role of 
the attractive forces on both the statics and the dynam- 
ics essentially becomes negligible and could be treated 
perturbatively. 

The important result obtained in our study is that the 
dynamics in region (II) at temperatures characteristic of 
viscous liquid behavior is strongly influenced, and in a 
highly nontrivial and nonperturvative way, by the attrac- 
tive forces. It seems therefore unlikely that an extrapola- 
tion from either regime (I) or regime (III) could fruitfully 
describe real glass- forming liquids in region (II). Thus, 
the common assumption that the behaviour of dense 
nonassociated liquids is determined by the short-ranged 
repulsive part of the intermolecular potentials dramati- 
cally breaks down for the dynamics in the viscous liq- 
uid regime. This striking result should motivate further 
work to reconsider the interplay between structure and 
dynamics, density and temperature, or else jamming and 
thermal activation in the glass formation. 
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